Scaling functions applied to three-body recombination of ^^^Cs 

atoms 

L. Platter^' and J. R. Shepard^'S 

^Department of Physics, The Ohio State University, Columbus, OH 43210 
^Department of Physics and Astronomy, 
Ohio University, Athens, OH 4 5701, USA 
^Department of Physics, University of Colorado, Boulder, CO 80309, USA 

(Dated: May 7, 2009) 

Abstract 

We demonstrate the implications of Efimov physics in the recently measured recombination rate 
of ^^'^Cs atoms. By employing previously calculated results for the energy dependence of the re- 
combination rate of '^He atoms, we obtain three independent scaling functions that are capable 
of describing the recombination rates over a large energy range for identical bosons with large 
scattering length. We benchmark these and previously obtained scaling functions by successfully 
comparing their predictions with full atom-dimer phase shift calculations with artificial ^He po- 
tentials yielding large scattering lengths. Exploiting universality, we finally use these functions 
to determine the 3-body recombination rate of ^'^^Cs atoms with large positive scattering length, 
compare our results to experimental data obtained by the Innsbruck group and find excellent 
agreement. 
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I. INTRODUCTION 



In atomic physics the term universahty refers to phenomena which are a result of a two- 
body scattering length a much larger than the range R of the underlying potential and 
do not depend on any further parameters describing the two-body interaction. The non- 
relativistic three-body system also exhibits universal properties if a ^ -R, but an additional 
three-body parameter is needed for the theoretical description of observables. Therefore, one 
three-body observable can be used (e.g. the minimum of the three-body recombination rate 
a=Ko) to predict all other low-energy observables of such systems. A particularly interesting 
signature of universality in the three-body system is a tower of infinitely many bound states 
{Efimov states) in the limit a = ±00 with an accumulation point at the scattering threshold 
and a geometric spectrum : 

E^") = [e-^^/^'^y^-^'h^l/m, (1) 

where is the binding wavenumber of the branch of Efimov states labeled by ri*. The 
three-body system displays therefore discrete scaling symmetry in the universal limit with 
a scaling factor factor e'^/*". In the case of identical bosons, Sq ~ 1.00624 and the discrete 
scaling factor is e'^^^° ~ 22.7. These results were first derived in the 1970's by Vitaly Efimov 
1, 2 1 and were rederived in the last decade in the framework of effective field theories (EFT) 

Recently, experimental evidence for Efimov physics was found by the Innsbruck group 
[sl. Using a magnetic field to control the scattering length via a Feshbach resonance, they 
measured the recombination rate of cold ^'^^Cs atoms and observed a resonant enhancement 
in the three-body recombination rate at a ~ — 850ao which occurs because an Efimov state 
is close to the 3-atom threshold for that value of a. The three-body recombination rate 
for atoms with large scattering length at non-zero temperature has been calculated with 
a numbej;_ of different models or based on the universality of atoms with large scattering 
lengths 



, [1^ . However, a striking way to demonstrate universality is to describe 



observables of one system with information which has been extracted from a completely 
different system. In jll|, the authors considered Efimov's radial laws which parameterize 
the three-atom S-matrix in terms of six real universal functions which depend only on a 
dimensionless scaling variable, x = {ma^E /h^Y^'^ , and phase factors which only contain 
the three-body parameter. In this work, simplifying assumptions justified over a restricted 



range of x were made to reduce the six universal functions required to parameterize the 
three-body recombination rate to just a single function. This function was then extracted 
from microscopic calculations of the recombination rates for ^He atoms by Suno et a/.[12|. 
In a recent paper, Shepard [13] calculated the recombination rates from atom-dimer elastic 
scattering phase shifts for four different ''He potentials (the so-called HDFB, TTY, LM2M2 
and HFDB3FCII potentials) and was able to obtain two universal functions. 



Here, we relax all but one of the simplifying assumptions made in 



ll| and extract a set of 



three independent universal functions capable of parameterizing the three-body recombina- 
tion rate over a wide range of energies. We test the performance of these universal functions 
using "data" generated from phase shift calculations [l3] employing artificial short-range ^He 
potentials. Finally, we use the new universal functions to calculate the scattering length and 
temperature dependent recombination rate for ^^^Cs atoms as measured by the Innsbruck 



group 



Sj and comment on our results. 



II. THREE-BODY RECOMBINATION 



Three-body recombination is a process in which three atoms collide to form a diatomic 
molecule (dimer). If the scattering length is positive and large compared to the range of 
the interaction, we have to differentiate between deep and shallow dimers. Shallow dimers 
have an approximate binding energy of -Eshaiiow — ^^/(ma^) <^ h^/{mR^). The binding 
energy of deep dimers cannot be expressed in terms of the effective range parameters and 
-E'deep ~ /{mB?). If the underlying interaction supports deep bound states, recombination 
processes can occur for either sign of a. In a cold thermal gas of atoms, recombination 
processes lead to a change in the number density of atoms ua 

jnA = -L^n\, (2) 

where L3 denotes the loss rate constant. The recombination coefficient, to which L3 is 
proportional, can be decomposed into 

K^{E) = irshallow(^) + KdeepiE) , (3) 

and the recombination rate into the shallow dimer can be further decomposed into contri- 
butions from the channels in which the the total orbital angular momentum of the three 



3 



atoms has a definite quantum number J according to 

oo 

i^shaUow(^)=5^i^^^n^) • (4) 
J=0 

For now, let us consider recombination via the shallow dimer only. If the collision energy E 
is small compared to the natural energy scale / {mB?), the recombination rate fTshaiiowl-E) is 
a universal function of the collision energy E , scattering length a and three-body parameter 
a^Q. The universal function depends on the dimensionless scaling variable defined as 

X = {ma^E/nY^ . (5) 

For J > the recombination rate does not depend on the three-body parameter a^^o and the 
implications of universality are therefore particularly simple, namely 

K^-^) = fj{x)haym . (6) 

However, K^^^ depends log-periodically on a*o (this is the signature of Efimov physics!) and 
is related to the S-matrix for elastic atom-dimer scattering through 

K^'\E) = - \Sad,ad\') , (7) 

Efimov's radial law then gives the dependence on complex universal functions and the three- 
body parameter a^^o which defines the scattering length for which the recombination rate 
has a minimum as 

„2 /™\ „2jso ln(a/a*o) 
o / \ 12 V / / o\ 

SaDAD = S22{X) + i_3^^(^),2..oln(aMo) " 

The functions Sn and Si2 are known at threshold 



-27rso 

5 



sii(O) = -e 
si2(0) = Vl - e-4^«oe^''°° , 

S22(0) = e2^^°°e-2"^° , (9) 

with 6oo = 1.737. It follows that |sii(0)| ~ 0.002. The first simplifying assumptions being 
made in was that this function remains small {i.e.; <^ 1) for all x and can be ignored. 
Then the energy dependent recombination rate can be written as 



K^"\E) = l^^^\l-{rl,-rt, + 2r,,rl,cos[<^ + 2so\og{a/a,o)]) 
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ha 
J m 

(10) 



where we have set Sij = Vij exp{i(j)ij) and $ = 022 — 20i2. Under the assumption that sn 
can be neglected the recombination rate depends therefore on the three real- valued function 
ri2(x), r22(a;) and It is worth noting that the expression in Eq. fllOl) is symmetric 

under exchange of ^-^d r22. However, the threshold conditions in Eq. can be used to 
to attribute the correct fit solutions to the universal function. 

As also discussed in Ref . jl^ , the effects of deep dimers can easily be incorporated through 
one additional parameter rj^ by making the substitution 

In a^o In a^o - iV*/-^o (H) 
unitarity the resulting effect on the recombination into shallow 



in, e.g. Eq. (|T0|) . Employing 



dimers can be written as 



li 



K^lnoJE) = ^^^(l - + .,2(x)^e— - (1 - e--)|.,2(x)r) ^ . 

(12) 

Note that in deriving this expression we assumed again that Su ~ 0. In the same manner 
one can derive an expression for the recombination rate into deep dimers 

K,ee,{E) = ^^^^(l-e-^''*)(l-|si2(a:)p) — . (13) 
III. ALTERNATIVE PARAMETERIZATIONS 



Starting with S-matrix element for 3-atom to dimer-atom scattering, it was shown in 11 1 
that under the assumption su = the recombination rate can be written as 

a 



sin[so ln( — )] (l + hi{x) + ih^{x)) 

+ cos[soln(— )](/i2(a;) +i/i4(x)) 
0*0 



(14) 



^ ■)\(I12(X)+Ihi(x)] ) 

m 

where Cmax ~ 67.1 and the hi are real- valued functions of x. Additionally,it was assumed 
that the imaginary part of the above amplitude can be neglected 

sin[soln(a/a^,o)](l + /ii(a;)) + cos[soln(a/a*o)]/i2(a;)| ha^/m. (15) 

This is well justified by direct calculations of the J = recombination rates for ^He atoms 
which display pronounced minima at approximately -^breakup — 20 mK and which can 




FIG. 1: The scaling functions hi (solid line) and /12 as a function of x. 



be explained by this assumption. Then the functions /13 and /14 can be set to in Eq. (ITUI) . 
The resulting expressions were employed in [3] to extract hi and /12 for x < 1.1. Although 
hi and h2 were determined by fitting to values of K^^\E) calculated using just two of the 
four atom-atom potentials considered, they were found to accurately account for the results 
for all 4 potentials as expected from universality. We have recalculated the /i-functions using 
the results for the three-body recombination obtained using the LM2M2 and HFDB3FC1I 
potentials and have fitted a polynomial to our results over the energy range < a; < 1.2 

hi{x) = -0.0234437X + 0.0550298a;2 - 1.03776x^ + 1.18985x^ - 0.471592x^ , 

h2{x) = 0.0338266a; - 0.233836^^ + 0.182564a;^ - 0.0895055a;^ + 0.0461793a;^ . (16) 

The functions are displayed in Fig. [H The effect of deep dimers on the recombination 
rate into the shallow dimer can easily be incorporated by making the substitution In a*o ~^ 
lna*o — iv*/-^o in Eq. ([Til) 



K^^\E) = Cmax cosh^r7*(sin[soln(a/a*o)](l + /ii(a;)) + cos[soln(a/a*o)]^2(a^))' 



-Fsinh^r7*(cos[soln(a/a*o)](l + hi{x)) - sin[so ln(a/a*o)]^2(a;)) -(17) 



ha'^ 



m 



To take the effects of the recombination rate into deep dimers into account it was assumed 



m 



ll| that Kdccp{E) is a function varying slowly with energy and that it can therefore be 
approximated with 



C 



deep 



(1 



ha 



m 



(18) 
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FIG. 2: The exact recombination rates and the corresponding results obtained with scahng (sohd 
hnes) and universal functions (dashed lines) of the HFDB (circles) and TTY (triangles) potentials. 

IV. EXTRACTION OF THE UNIVERSAL FUNCTIONS 

By fitting Eq. fnU]) to the recombination rates of all four ^He potentials, we were able 
to determine the functions ri2(x), r22(x) and ^{x). Our results are smooth functions for 
X > 0.2 and the radial functions approach the known threshold values from Eq. (Q for 
decreasing x. For x < 0.2, we are not able to find a reliable fit which is indicated by the 
rapid variation of the function $ in Fig. [3] in this region. 

To display the qualities of our fit we compare the exact recombination rates obtained with 
the TTY and HDFB potentials to the rates calculated with the newly obtained universal 
functions These results for these functions are displayed in Fig. [2l This figure contains 
also the recombination rate obtained with the /i-functions. While the new set of universal 
functions seem to provide slightly better results for the HFDB potential at larger energies, 
the /i-functions perform equally well for these potentials at lower energies. 

To test our new parameterizations we have generated three artificial potentials (which 
we call I, II and III) characterized by different three-body parameters a^^Q (with a/a^^Q = 
1.384,1.188 and 1.780, respectively) but having approximately the same ratio of R/a as 
the real ^He potentials used in this work. We have calculated the recombination rates for 
these potentials and use these results to benchmark our universal functions. Our results are 
displayed in FigHl We find that the new set of functions is capable of describing the recom- 
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FIG. 4: The exact recombination rates and the corresponding results obtained with the /i- functions 
(solid lines) and the universal functions (dashed lines) for the potentials I (circles), II (squares) 
and III (triangles) potentials. 

bination rates of these potentials over a relatively large range of x. Again we benchmark 
also the rates obtained with hi and h2 and find that this set of scaling functions describes 
the exact results better at x < 1 than the scaling functions ri2(x), r22{x) and $(x). This 
is surprising at first sight since one certainly expects to obtain a better description of the 
recombination rate with three instead of two functions. We speculate that the functional 
form in Eq. ffTOj) results in stronger constrains on the universal functions than Eq. f|T^ does 
on the h functions. All the potentials, however, contain finite range effects which are not 
accounted for in Eq. (IHl). It is therefore very likely that better fits - using the same ap- 
proximation - can be obtained from recombination rates calculated in the exact zero-range 
limit. 



V. RESULTS FOR CESIUM 



In the previous section we found that we can obtain a very good overall description of the 
recombination rate of systems with a large scattering length if we employ the functions hi 
and h2 for energies smaller than -Eshaiiow and the universal function ri2{x), r22{x) and $(a;) 
for energies larger than -Eshaiiow Using these functions at energies close to the minimum in 
the recombination guarantees a more appropriate treatment of the effect of deep dimers on 
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FIG. 5: The 3-body recombination length for Cs for a=KO = 210ao and three different values 
of the parameter ry^KiO (solid line), 0.01 (dashed line), and 0.06 (dotted lines) plotted together with 
the experimental results of the Innsbruck experiment (triangles) 



the recombination rate, which are expected to have the largest effect in this region. 

The form of the functions fj{x) in Eq.(l6]) and therefore the contribution to the recombi- 
nation from channels with higher total angular momentum J has been previously analyzed 



m 
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13| , we thus take these channels into account by using appropriate parameterizations 



for the functions fj{x). ^^^Cs atoms can recombine into deep and shallow dimers. As men- 
tioned above, a deep dimer is so strongly bound that it cannot be described within the EFT 
for short-range interactions as the binding energy is larger than /{mR^). We account for 
such processes by letting lna*o lna*o — iv*/^o as also discussed above. We then calculate 
the temperature dependent recombination rate by calculating 

dE e-'^/^'''^^^ KsjE) 
- Qjl^dEE^e-^/C'BT) ■ (19) 

The weight factor E"^ comes from using hyperspherical variables for the Jacobi momenta. 

/ \ 1/4 

In Fig. [5] we display our results for the recombination length = f^^j of ^^^Cs 
atoms. It can be seen that the results agree very well with the experimental results obtained 
by the Innsbruck group at T = 200 nK. 



10 



VI. SUMMARY 



In this paper we have used the results from different ^He atom-atom potentials to extract 
and to test the predictive power of universal functions. In doing so, we have relaxed all 
but one simplifying assumptions which was made in previous work P, Q. We have de- 
termined a third universal scaling function which allows for a description of the three-body 
recombination rate of systems with large scattering length over a greater range of breakup 
energies. 

We have tested the quality of our parameterizations with artificial finite range potentials 
which are appreciably different from the original Helium potentials but which display uni- 
versal effects in three-body sector. We have found that our three real universal functions can 
describe the recombination of these artificial potentials reasonably well which gives further 
evidence that the assumptions made in ll| were well justified. We also found, however, 
that the previously calculated scaling functions hi and h2 give an overall better description 
of the recombination rate for energies E < -Eshaiiow The scaling functions hi and h2 which 
can be represented analytically with a simple polynomial fit given in Eq. flTB]) are therefore 
a useful tool to test recombination rate calculations for systems with large scattering length. 

Finally, we have used both sets of universal functions together to compute the recombi- 
nation length for ^^'^Cs atoms for different values of the parameter r^^, which approximately 
accounts for the effect of deep dimer states and have compared our results with experimental 
data obtained by the Innsbruck group[5|]. 

Although our results show very good agreement with the data, sensitivity to ?7=k is insuf- 
ficient to permit a precise determination of this parameter. Overall, we consider our results 
to be an excellent example of how few-body systems with large scattering length exhibit uni- 
versal features. The low-energy properties of ^He atoms allow us to compute accurately the 
low-energy properties of a gas of a completely different element, ^^^Cs, which at first glance 
has little in common with ^He. Nevertheless, we point out that the results cannot be thought 
of as complete treatment of the problem at hand. For example, not only did we make the as- 
sumption that Sii does not contribute significantly to the recombination coefficients, we also 
extracted the functions from data sets obtained with finite range potentials. Although the 
impact of range corrections is known to be small for realistic Helium atom-atom potentials 
as -R/a ~ 0.1, it needs to be pointed out that range corrections are expected to be sizable 
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for large enough energies. To obtain all universal functions Sij relevant to the recombination 
rate, a calculation in the limit i? — > seems therefore to be necessary ^. Furthermore, it 



is already understood how to include range corrections systematically in the framework of 

14, [is], Q- Indeed, this approach has already been used to calculate 

n " 



effective field theory 



range corrections to the recombination rate into a shallow dimer Il9l |. Thus, further 
effort should be devoted to include these effects in the calculation of the energy- dependent 
recombination rate. 
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